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Let  a  moving-average  interpolation  formula  for  equally  spaced  data,  exact 
for  the  degree  r,  have  a  basic  function  L  e  C®’1  of  finite  support  with 


•  (») 

a 

/V 


piecewise  continuous*  Such  a  formula  is  called  "smoothest"  when  the 


integral  of  the  square  of  over  the  support  of  L  is  smallest.  If  m, 

a  » 

r,  and  the  support  of  L  are  given,  either  there  is  no  such  formula  or  there 
is  a  unique  smoothest  formula,  for  which  L  is  a  piecewise  polynomial  of 
degree  at  least  r  and  at  most  max(r,  an  -  1),  uniquely  characterised  by 
certain  conditions  on  the  location  of  its  knots  and  the  jumps  occurring 
there .  A  similar  result  is  obtained  if  consideration  is  limited  to  formulas 
that  preserve  (i.e.,  do  not  smooth)  the  given  data. 
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SIGNIFICANCE  AND  EXPLANATION 


In  local  moving-average  interpolation  of  equally  spaced  data,  each 
interpolated  value  is  calculated  as  a  weighted  average  of  a  few  given 
■>r  inates  situated  near  the  ordinate  that  is  being  approximated.  The  weight 
applied  depends  only  on  the  distance  between  the  argument  of  the  ordinate  to 
which  the  weight  is  applied  and  the  argument  of  the  interpolated  value.  An 
example  is  a  procedure  that  was  often  used  in  the  construction  of  mathematical 
tables  before  computers  became  available.  Every  fifth  or  tenth  value  was 
calculated  from  a  series  expansion,  and  the  intermediate  values  were  obtained 
by  local  Newton-La grange  interpolation.  This  amounts  to  fitting  a  piecewise 
polynomial  function.  The  piecewise  curve  so  derived  had  "corners,"  but  the 
discontinuities  in  the  first  derivative  were  too  small  to  be  of  any  importance 
when  a  smooth  mathematical  function  was  being  interpolated. 

In  the  nineteenth  century  British  actuaries  noted  that  the  "corners"  in 
the  curve  of  interpolated  values  were  often  objectionable  when  empirical  data 
were  being  interpolated,  and  they  published  a  number  of  local  interpolation 
formulas  in  which  piecewise  polynomials  of  higher  degree  were  fitted,  and  the 
additional  degrees  of  freedom  so  obtained  were  utilized  to  secure  smooth 
junction  of  adjoining  polynomial  arcs.  In  1946  Schoenberg  showed  that  a  local 
moving-average  interpolation  formula  is  fully  characterized  by  a  certain 
function  of  finite  support  having  a  bell-shaped  graph,  similar  in  appearance 
to  a  probability  distribution,  except  that  it  usually  assumes  some  negative 
values  in  the  tails.  This  he  called  the  basic  function  of  the  interpolation 
formula .  In  1 954  the  present  authors  published  a  paper  in  which  Schwartz 
distributions  and  Schoenberg's  basic-function  concept  were  used  to  develop  a 
general  theory  of  smooth-junction  local  interpolation  formulas. 

If  we  fix  the  support  of  the  basic  function,  the  degree  (of  polynomials) 
for  which  the  formula  is  to  be  exact,  and  the  order  of  derivatives  to  be  used 
in  judging  smoothness,  it  is  shown  in  the  present  paper  that  there  is  then  a 
unique  interpolation  formula  of  the  class  so  defined  that  is,  in  a  certain 
sense,  smoothest.  Two  cases  are  considered:  that  in  which  the  curve  of 
interpolated  values  is  required  to  pass  exactly  through  the  given  data  points, 
and  the  more  general  case  in  which  greater  smoothness  is  obtained  by  dropping 
i:his  requirement. 
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SMOOTHEST  LOCAL  INTERPOLATION  FORMULAS  FOR  EQUALLY  SPACED  DATA* 

T.  N.  E.  Greville  and  Hubert  Vaughan1 

1 «  INTRODUCTION .  Schoenberg  pointed  out  in  1946  (7)  that  a  large  class 
of  local  interpolation  formulas  for  equally  spaced  data  can  be  expressed  In 
the  form 


v„  -  [  L(x  -  v)y  ,  (1.1) 

V—  V 

where  yv  denotes  a  given  ordinate#  vx  is  an  interpolated  value#  and  L(x) 
is  a  given  function  called  by  Schoenberg  the  basic  function  of  the  interpola¬ 
tion  formula.  This  class  includes  the  numerous  formulas  of  so-called  "oacula- 
tory  interpolation"  published  by  actuarial  writers  (for  additional  references 
see  [3]).  For  the  latter  formulas  L  is  typically  a  piecewise  polynomial 

1  O 

function  of  finite  support  belonging  to  continuity  class  C  or  C  . 

Also  included  is  what  may  be  called  moving  Newton-La grange  interpolation# 
often  used#  before  aomputers  became  available,  in  the  preparation  of  tables  of 
mathematical  functions.  An  example  would  be  the  case  in  which  the  function 
f  is  interpolated  in  (Vh,  (v  +  1)h),  v  being  any  integer,  by  means  of  the 
cubic  p  uniquely  determined  by  the  four  oonditions 

p(x)  -  f(x)  (x  -  (v  +  j)h»  j  -  -1,0, 1,2)  . 


* 

A  6-page  synopsis  of  results,  without  proofs,  appeared  in  Approximation 
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In  this  example  L  is  a  continuous  piscswiss  cubic  with  support  in  (-2*2) < 
whose  first  derivative  is  discontinuous  at  x  •  -2,-1, 0,1, 2  (see  Fig.  1). 
These  discontinuities  did  not  give  rise  to  any  problems  in  interpolating 
smooth  mathematical  functions,  but  are  undesirable  if  one  is  interpolating 
empirical  data. 

Formula  (1.1)  is  called  reproducing  when  L  is  such  that  vy  -  y  y  for 
every  integer  v,  whatever  may  be  the  values  of  the  quantities  yv.  It  is 
clear  that  (1.1)  is  reproducing  if  and  only  if 

L(\>)  *  80v  (V  *  •••#-1*0,1,...)  .  (1.2) 

where  is  a  Kronecker  symbol.  In  practice  an  interpolation  formula  that 

is  not  reproducing  smoothes  as  well  as  interpolates,  since  each  given  ordinate 
yv  is,  in  general,  replaced  by  an  adjusted  value  vy.  (Whether  the  adjust¬ 
ment  does,  in  fact,  actually  increase  the  smoothness  of  the  data  depends  on  a 
judicious  choice  of  Li  see  [7]  .) 

Figure  1  shows  the  graphs  of  three  typical  basic  functions.  Note  that 
Karup's  formula  and  the  Newton-La grange  central  third-difference  formula  are 
reproducing,  while  Jenkins'  "modified"  third-difference  formula  is  not.  On 
the  other  hand,  note  that  the  Newtonian  graph  has  oorners,  while  the  others  do 
not. 

Formula  (1.1)  is  called  exact  for  the  degree  r  when  L  is  such  that 
the  formula  gives  exact  values  whenever  it  is  used  to  interpolate  a  polynomial 
of  degree  r  or  less.  In  other  words,  using  to  denote  the  class  of 

polynomials  of  degree  r  or  less,  L  is  such  that,  for  every  p  6  i^, 
y  *  p(v)  for  all  integers  v  implies  v  »  p(x)  for  all  real  x. 

V  X 


Figure  1 

In  the  case  of  moving  Newton-La  grange  interpolation,  r  is  merely  the 
degree  of  the  polynomial  arcs  employed.  For  the  actuarial  formulas,  r  is 
less  than  the  degree  of  the  piecewise  polynomial  function  L,  the  "degrees  of 
freedom"  thus  gained  being  utilized  to  increase  the  order  of  continuity  of 
L.  The  latter  is,  of  course,  also  the  order  of  continuity  of  the  composite 
interpolating  function,  for  if  L  ©  C*"”1,  then  it  is  clear  from  (1.1)  that 
v  e  c”*-1. 

When  L  is  discontinuous  (as  occurs,  for  example,  in  the  case  of 
symmetrical  moving  Newton-Lagrange  interpolation  of  even  degree),  the 
definition  of  exactness  for  the  degree  r  requires  interpretation.  In  such  a 


case,  it  must  be  assumed  that  L,  though  discontinuous,  is  nevertheless  such 


that,  for  avary  p  ®  «r# 

m 

d(x)  -  I  L(x  -  w)p(v) 

has  only  ranovabla  dlsoontlnultlas,  and  thay  ara  removed  by  taking 

rf(x)  -  #(x  ♦  0)  -  #(x  -  0)  . 

2 .  HINIHIZED-DIRIVATIVI  fORMPLAB .  Lat  I  -  (<*,B)  ba  a  finita  opsn 
Interval  on  the  real  line,  and  let  FIrm  denote  the  eet  of  interpolation 
formulas  of  the  form  (1.1)  that  are  exact  for  the  degree  r,  and  have  a  basic 
function  I.  6  C*-1  with  its  support  contained  in  x  and  with  t/** 
piecewise  continuous •  Also  let  FjJJ  denote  the  subset  of  rirm  consisting 
of  reproducing  formulas.  It  follows  from  (1.2)  that  !•  •Pty  unless 

0  e  X.  By  a  piecewise  continuous  function  we  mean  one  having  only  jump 
discontinuities  and  at  most  a  finite  set  of  these. 

In  each  of  the  classes  FInB  or  FjJjj  **  find  that 

formula  which  is  in  some  sense  smoothest.  MS  shall  judge  smoothness  by  the 
closeness  to  zero  of  the  mth  derivative  of  the  interpolating  function  vx» 

Mow,  m-fold  differentiation  of  (1.1)  gives 

m 

y(m)  m  l  L(m)(x  _  v)  (2.1) 

v— • 

almost  everywhere .  As  we  have  some  latitude  in  the  choice  of  the  basic 
function  L,  but  none  as  regards  the  given  ordinates  yy,  (2.1)  suggests 
that  the  values  of  vxra^  will  be  closer  to  zero  than  would  otherwise  be  the 
case,  if  L  is  chosen  so  that  the  values  of  L are,  in  some  sense,  as 
close  to  zero  as  possible.  Accordingly,  we  shall  call  a  given  formula  of  one 


of  th#  class**  fjrB  or  *  »lnlal»*d-d*rlv*tlva  formula  (»df)  of  it* 

class  if  th*  quantity 


In  other  word*,  we  have  shown  that,  tor  a  given  M,  by  minimising  J  we 
minimise  «  certain  upper  bound  to  the  integral  over  a  unit  interval  of  the 
abeolute  value  of  v£*)  • 

Minis  iced- derivative  interpolation  formulae  were  previously  defined  by  us 
in  (3],  and  a  few  example*  were  given,  but  no  general  theory  we*  developed. 

In  this  paper  we  shall  show  that  the  class  fIrv  1*  empty  for 
u  <  r  ♦  1 ,  i*  empty  or  contain*  a  single  formula  (of  moving  Me wton-Le grange 
interpolation)  when  V  •  x  ♦  1,  and  i*  infinite  «dten  u  >  r  ♦  1.  In  the 
latter  case,  w*  shall  show  that  there  i *  a  unique  mdf,  and  shall  characterise 
this  formula  in  a  way  that  lead*  to  an  algorithm  for  it*  determination  in  any 
particular  case. 

It  was  previously  noted  that  the  class  is  empty  unless  the  open 

interval  I  contains  the  origin,  then  this  condition  is  satisfied,  w*  shall 
find  that  is  identical  to  PIn|  when  the  number  of  Integers  contained 

in  I  does  not  exceed  r  ♦  1.  This  implies,  of  course,  that 
r  ♦  1  <  W  <  r  ♦  2,  but  the  aonverse  is  not  true. 

when  the  number  of  integers  contained  In  I  exceeds  r  ♦  1  (and  0  is 
among  them),  F^rm  ls  *  proper  subset  of  Fjm,  and  the  former  contains  a 
unique  mdf  different  from  that  associated  with  the  latter.  In  this 
(reproducing)  case  too,  we  shall  obtain  a  characterisation  of  the  mdf  leading 
to  an  algorithm  for  its  determination. 

In  order  to  arrive  at  the  results  just  described,  it  is  first  necessary 
to  express  the  requirement  that  (1.1)  be  exact  for  the  degree  r  in 
manageable  form  as  a  set  of  constraints  on  the  basic  function  L.  A 
digression  for  this  purpose  is  the  subject  of  the  next  section. 
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3.  MAINTCTANCg  OF  DBGREt .  Schoenberg  noted  in  (7)  thet  the  implications 
of  exactness  of  formula  (1.1)  for  the  degree  r  beooes  clearer  if  ooneldered 
in  relation  to  a  certain  weaker  oonditlon.  This  weaker  aondition.  in  a 
modified  form,  was  utilised  by  us  in  (3)  and  is  used  again  here,  let  l  be  a 
given  function  with  its  support  contained  in  I.  and  let  a  function  p  and  a 
function  d  be  related  by  the  formula 

m 

d(x)  •  l  H(x  -  v)p(u)  .  (3.1) 

This  relation  may  be  regarded  as  a  transformation  TH  that  transforms  p 
into  d,  or 

d  ■  T||P  . 

It  is  evident  that  TH  is  a  linear  operator,  he  shall  say  that  T„ 
maintains  the  degree  r  if  it  maps  the  specs  into  Itself,  or.  in  other 

words,  if  p  @  implies  d  9  *r  • 

An  important  special  case  of  maintenance  of  degree  is  that  in  tdtich  TM 
annihilates  *  :  in  other  words,  d  is  ldsntlcally  zero  whenever  p  9  *r . 

Schoenberg  in  [7]  defined  a  transformation  TH  that  prsssrves  the 
degree  r  as  one  having  the  property  that,  for  every  p  •  *f,  d  le  e 
polynomial  strictly  of  the  same  degree  as  p  with  the  same  leading 
coefficient.  He  showed  that  TH  has  this  property  if  ite  characteristic 
function  (Fourier  transform  of  the  basic  function  H)  has  the  value  1  for  the 
argument  0  and  zeros  of  order  r  ♦  1  for  all  nonvanishing  integral  multiples 
of  2*  • 

We  showed  in  (31  that  if  H  has  an  (r  ♦  1)th  derivative  in  the  sense 
of  distributions,  then  TH  maintains  the  degree  r  if  and  only  if  the 
convolution  of  that  (r  ♦  1)th  derivative  with  every  element  of 
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vanishes.  A  more  limited  result ,  that  can  be  stated  and  proved  without 
introducing  distributions ,  will  auffioe  here,  and  is  contained  in  the  newt 
theorem . 

If  H  has  a  piecewise  continuous  jth  derivative,  let  b^  denote  the 
jueip  of  at  x  -  t. 

Theorem  3.1 .  Xf  Ty  main tain a  the  degree  r  and  if  M  is  pieoewise 
continuous  and  has  piecewise  continuous  derivatives  of  all  orders,  then,  for 
all  real  t. 


l  vlh 

i*1”' 


(1  •  0,1,  ...,r»  J  •  0,1,...) 


0.2) 


If  H  is  a  piecewise  polynomial  function  with  finite  support  and  satisfies 
(3.2),  then  T„  maintains  the  degree  r. 

Proof.  If  ws  take  p(v)  •  v*  in  0*1),  the  left  member  of  0.2)  is  the 
jump  of  d^J(x)  at  x  •  t.  But,  if  TH  maintains  the  degree  r  and  i  is 
one  of  the  Integers  0,1,..., r,  4  is  a  polynomial  and  this  jump  vanishes. 

Thus  (3.2)  is  established. 

On  the  other  hand,  let  H  be  a  piecewise  polynomial  function  of  finite 

support  satisfying  (3.2),  and  let  p  9  i  in  0.1).  Then  4  and  all  its 

derivatives  are  continuous  everywhere.  But,  if  d  is  the  maximum  degree  of 

the  polynomial  arcs  composing  H,  then  vanishes  almost  everywhere  by 

(3.1).  Since  is  continuous  everywhere,  it  is  therefore  identically 

zero.  It  follows  that  4  6  *.• 

a 

If  S  denotes  the  "central-difference"  operator  defined  by 


st <x)  -  f<x  ♦  y2  )  "  f (x  -  V2  )  » 


we  have  a leo  f  rom  (3.1) 


I  f 


8r>1#(x)  •  I  p<v)5r^1H(x  -  v) 


Expanding  5r+1H(x  -  v)  In  terms  of  H(x)  values  by  tha  wall  known  binomial 
formula  and  rearranging  tarma  gives 

«r+1#(x)  -  I^Htx  -  v)«r+,p(v>  . 

whara  tha  luwitlon  is  ovar  all  tha  Integers  whan  r  ia  odd,  and  over 

all  tha  real  numbers  of  tha  form  integer  ♦  V2  whan  r  ia  even.  Note  that 
tha  required  rearrangement  of  terms  ia  permissible  because  the  support  of  H 
is  finite . 

Now,  since  p  0  6r+1p<v)  “  0  for  all  v.  Therefore  6r+1d(x)  ■  0 

for  all  real  x.  But,  a  polynomial  whose  (r  ♦  l)th  derivative  vanishes 
identically  belongs  to  Therefore  TH  maintains  the  degree  r.  □ 

If  Th  maintains  tha  degree  r,  then  there  is  a  differential  operator 
of  order  not  exceeding  r,  which  we  shall  call  the  signature  of  TH  and 
shall  denote  by  S^,  of  the  form 

S„  -  1  a^D*  ,  (3.3) 

i-0 

that  is  equivalent  to  over  .  In  other  words,  S^p  ■  THp  whenever 

p  6  In  (3.3)  D  denotes  differentiation.  The  following  theorem  is  an 

immediate  consequence  of  the  preceding  definitions. 

Theorem  3.2.  TH  is  exact  for  the  degree  r  if  and  only  if  it  maintains 
the  degree  r  and  its  signature  is  the  identity  operator . 

Vife  can  express  SH  in  terms  of  H  in  various  ways.  Thus,  the 
coefficients  a^  of  (3.3)  are  given  by 


I  VAH(V)  (i  -  0,1,...,r>  .  (3.4) 

1  V«-°° 

If  H  is  integrable,  we  have  also 

i  ” 

/  x^fxjdx  (i  *  0,1, ...,r)  . 


ai  "  60i 


If  a  basic  function  L  satisfies  (1.2),  then,  by  (3.4),  the 
coefficients  a^  in  the  expression  (3.3)  for  SL  are  given  by 

(i  “  0,1, ...,r)  . 

Thus  we  have  established  the  following  corollary,  previously  noted  in  [3,7]  . 

Corollary  3.3.  The  interpolation  formula  (1.1)  is  exact  for  the  degree 
r  if  it  maintains  the  degree  r  and  is  also  reproducing. 

The  existence  and  properties  of  the  signature  SH  were  established  in 
[3]  (though  the  term  "signature"  does  not  appear  there)  using  the  concept  of 
distributions.  However,  what  has  been  stated  here  is  easily  verified  by 
elementary  means.  A  similar  remark  applies  to  the  following  lemma,  in  which 
we  take 


p  =  .I(r  +  i)  . 

The  lemma  can  be  verified  by  noting  (after  some  algebraic  manipulation)  that 
the  function 


G(x)  *  l  (rtj)  K(x  -  P  -  j)  (3.5) 

j-0  3 

has  the  required  properties . 

Lemma  3 .4  .  If  K  is  piecewise  continuous  with  its  support  contained 

in  I  and  T-  annihilates  x  ,  then  there  exists  a  piecewise  continuous 
^  r 
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(3.6) 


function  G  with  support  contained  in  (#  +  p,  S  •  p),  such  that 

fir+1G(x)  -  K(x) 

for  all  real  x  such  that  the  left  member  is  defined. 

4.  CHARACTERIZATION  OF  GENERAL  mdf’s.  We  note  that  a  formula  (1.1)  that 
is  exact  for  the  degree  r  must  satisfy 

00 

l  ^Ux  -  v)  -  x1  (i  -  0,1,  ...,r)  .  (4.1) 

v=»-® 

If  the  support  of  L  is  aontained  in  I,  all  but  a  finite  number  of  the 
coefficients  L(x  -  v)  vanish  automatically.  If  W  <  r  +  1,  it  follows  that 
there  is  some  interval  for  x  within  which  each  of  the  r  +  1  linearly 

2  T* 

independent  functions  1,x,x  ,...,x  is  expressible  as  a  linear  combination 
of  r  or  less  given  functions.  This  is  impossible.  Therefore  Flrm  is 
empty  for  P  <  r  +  1 . 

If  p  =  r  +  1 ,  then  for  every  x  such  that  x  -  a  and  x  -  0  are 

nonintegers,  (4.1)  can  be  regarded  as  a  system  of  r  +  1  linear  equations  in 

the  r  +  1  unknown  values  of  L(x  -  V).  Moreover,  the  determinant  of  the 
matrix  of  coefficients  of  the  linear  system  is  a  Vandermonde,  and  therefore 
nonvanishing.  Thus,  the  system  has  a  unique  solution.  Now,  it  is  evident 
that  the  equations  are  satisfied  by  the  r  +  1  fundamental  functions  of 
Lagrange  interpolation  (or  extrapolation)  for  the  function  value  corresponding 
ro  thd  argument  x,  given  those  corresponding  to  the  r  +  1  arguments  v 
for  which  L(x  -  v)  is  undetermined.  Moreover,  each  of  these  fundamental 
functions  is,  indeed,  a  function  of  x  -  v,  as  (4.1)  requires. 

In  this  case  of  V  =  r  +  1 ,  L  is  discontinuous  at  those  arguments  x  e  I 

that  differ  by  an  integer  from  a  or  0,  except  in  the  special  case  in  which 
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a  and  8  are  themselves  Integers  and  also  061*  Only  In  this  special 
case  is  the  formula  reproducing  and  L  continuous  everywhere . 

we  conclude  from  the  preceding  discussion  that  the  class  FIrm  is  empty 
for  p  <  r  +  1 ,  and  also  for  P  ■  r  +  1  and  m  >  1 ,  while  for  P  *  r  +  1 
and  m  -  0  or  1 ,  it  is  either  empty  or  contains  a  single  formula . 

For  P  >  r  +  1,  *Irn  oontains  an  infinite  number  of  formulas  for  every 
nonnegative  number  m,  and  among  them,  as  we  shall  see,  a  unique  mdf .  The 
following  theorem  is  the  key  to  the  characterization  of  this  unique  mdf. 

Theorem  4.1.  For  any  nonnegative  integers  r  and  m,  and  for 
p  >  r  +  1 ,  the  class  FJrB  oontains  a  single  formula  whose  basic  function 
L  satisfies  the  following  three  conditions : 

(i)  L  is  a  piecewise  polynomial  function  of  degree  at  least  r  and  at 
most  d  =  max(r,2m  -  1). 

(ii)  Each  knot  of  I<  is  an  argument  that  differs  by  an  integer  from  a 
or  B  (or  both)  . 

(iii)  The  piecewise  polynomial  function  $  L  is  given  in 

(a  +  p,  (J  -  p)  by  a  simple  polynomial  of  degree  at  most  2m  -  1 . 

This  theorem  requires  interpretation  for  m  «  0.  In  that  case,  we 
interpret  a  polynomial  of  degree  -1  (in  condition  (iii))  to  mean  one  that  is 
identically  zero . 

We  shall  postpone  the  proof  of  this  theorem,  as  it  will  become  easier 
after  we  have  developed  some  further  paraphernalia.  However,  without  waiting 
to  prove  it,  we  shall  proceed  to  demonstrate  its  connection  with  the  existence 
of  a  unique  mdf.  For  this  purpose  we  shall  need  the  following  leirnna. 

Lemma  4 .2 .  Let  K  @  d"”1,  with  piecewise  continuous,  have  its 

support  in  I,  and  let  TR  annihilate  *  •  Let  H  be  piecewise  continuous, 
with  piecewise  continuous  derivatives  of  orders  1  to  m,  let  T^  maintain 
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IT  ^  1 

the  degree  r,  and  let  8  H  be  given  in  (a  +  p,  3  -  P)  by  a  simple 
polynomial  of  degree  at  most  2m  -  1 .  Then, 


/  H(m) (x)K(m) (x)dx  -  0  .  (4.2) 

Proof .  By  Lemma  3.4,  there  exists  a  function  G  0  d®”1,  with  support  in 
(a  +  p,  3  -  p),  such  that  is  piecewise  continuous  and  (3.6)  holds. 

Denoting  by  o  the  left  member  of  (4.2),  we  have 

O  -  /"  . 

If  8r+1G^m* (x)  is  expanded  in  terms  of  G*m*(x)  values,  the  finite  support 
of  G*m'  then  permits  rearrangement  of  terms,  so  that 

am  (-l)r+1  /  G(m)(x)6r+1H(m)(x)dx  . 

Since  G*m*  vanishes  outside  of  (a  +  p,  3  -  p),  and  6r+1H*m*  is  given  in 

that  interval  by  a  polynomial  of  x  ,  say  q,  we  have 

in—  1 

3-p 

0=  (-1)r+1  /  G(m) (x)q(x)dx  .  (4.3) 

a+p 

As  G  e  d®"1,  m-fold  integration  by  parts  now  gives  a  -  0,  as  required.  □ 
Theorem  4.3.  The  unique  interpolation  formula  determined  by  Theorem  4.1 
is  the  unique  mdf  of  the  class  FIrm» 

Proof  «  Let  L  be  the  basic  function  of  the  unique  formula  determined  by 
Theorem  4.1,  and  L^  the  basic  function  of  any  formula  of  FIrm.  Also  let 
J  and  J1  denote  the  corresponding  values  of  the  quantity  given  by  (2,2), 
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and  let  be  defined  by 


L,(x)  -  L(x)  +  (x)  .  (4.4) 

Then,  it  is  easily  verified  that  fulfills  the  conditions  required  of  K 

in  Lemma  4.2.  Similarly,  L  fulfills  the  requirements  for  H  in  that 
lemma .  Therefore ,  by  Lenna  4.2, 

00 

/  L(m) (x)kJb) <x)dx  -  0  .  (4.5) 


From  (4.4)  and  (4.5)  we  have 

00 

J1  -  J  +  /  [K<m) (x)]2dx  . 

_O0 

It  follows  that  J  £  Jj.  Moreover,  equality  holds  only  if  icjm^  vanishes 
almost  everywhere.  Therefore,  in  this  case,  is  a  step  function. 

But,  since  e  C1"-1,  is  continuous,  and  therefore  is 

identically  zero.  It  follows  that  K.  e  w  ..  But  a  polynomial  with  finite 
support  is  identically  zero,  and  so  L|  >  L.  □ 

5.  CHARACTERIZATION  OF  REPRODUCING  mdf's.  We  shall  first  dispose  of  the 
case  in  which  FIrm  and  Fj^P  are  identical.  Let  t  denote  the  number  of 
integers  contained  in  I. 

Theorem  5.1 .  If  FIrm  is  nonempty,  it  is  identical  to  Fj||P  if  and 
only  if  0  6  I  and  t  £  r  +  1 . 

Proof .  We  have  seen  that  FIrm  is  empty  for  v  <  r  +  1 .  Therefore  we 
must  have  w  £  r  +  1.  This  implies  that  t  >_  r .  In  fact,  t  ■  r  occurs  only 
when  M  »  r  +  1  and  a  and  8  are  integers.  By  (3.4)  and  Theorem  3.2, 


the  r  +  1  relations 


(5.1) 


I  vLuv)  -  6  (i  -  0,1,  ...,r) 

must  be  satisfied.  These  may  be  regarded  as  a  system  of  linear  equations  in 
the  quantities  L(V)  for  the  t  integers  v  contained  in  I.  If 
t  =  r  -i-  1,  the  matrix  of  coefficients  is  square.  (In  the  special  case  of 
t  =  r,  one  of  the  end  values,  a  or  8,  may  be  included.)  The  matrix  is 
also  nonsingular,  because  its  determinant  is  a  Vandermonde  and  therefore 
nonvanishing.  Thus,  the  linear  system  has  a  unique  solution.  However,  if 
0  e  I,  it  is  evident  that  the  values  given  by  (1.2)  are  a  solution,  and 
therefore  the  unique  solution .  Hence  the  formula  is  reproducing  and  FIrm 
and  are  identical . 

As  previously  pointed  out,  there  is  no  reproducing  formula  if  001. 

Now,  let  I  be  such  that  t  >  r  +  1,  and  let  f^  be  a  formula  of  FInn.  It 
follows  from  Theorem  4.1  that  such  a  formula  exists.  If  f1  is 
nonreproducing,  the  theorem  is  established.  Otherwise,  let  be  the  basic 

function  of  f1  and  t  the  largest  integer  in  I.  Then,  for  some  e  in 
(0,1/2),  [T-r-1-e,T+e]C  I.  Now,  consider  the  interpolation 
formula  f2  under  which,  for  x  in  (A  -  e,  A  +  e)  for  every  integer  A, 

00 

vx  =  I  L  (x  -  v)yy  +  k(x  -  A  +  e)m(x  -  A  -  e)mAr+1yx_T  ,  (5.2) 

V  being  arbitrary,  while,  for  all  other  values  of  x,  vx  is  given  by  the 
summation  term  only.  Here  A  is  the  usual  finite-difference  operator. 
Evidently  f2  belongs  to  FIm.  However,  (5.2)  gives 

vA  -  yx  +  k(-1)  c  A  yx_T  , 

and  the  formula  is  clearly  nonreproducing  for  k  *  0 .  □ 
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The  analogue  of  Theorem  4.1  for  the  reproducing  cbm  la  the  following 
theorem . 

Theorem  5.2.  For  any  nonnegative  Integer  r  and  positive  integer  a, 
and  any  finite  interval  1  «  (a, 6)  containing  0  and  such  that  V  >  r  ♦  1 , 
the  class  contains  a  single  formula  satisfying  the  following  three 

conditions: 

(i)  L  is  a  piecewise  polynomial  function  of  degree  at  least  r  and  at 
most  d  -  max(r,2m  -  1). 

(ii)  Each  knot  of  L  is  either  an  integer  or  an  argument  that  differs 
by  an  integer  from  a  or  0  (or  both) . 

(iii)  The  piecewise  polynomial  function  6r+1L  is  given  in 

(a+p,  &  -  p)  by  a  spline  function  of  degree  2m  -  1  with  simple  knots. 
The  knots  of  6  L  in  (a+p,  0-p)  are  at  the  integers  when  r  is  odd, 

and  at  the  arguments  of  the  form  integer  +  1/2  when  r  is  even. 

As  in  the  case  of  Theorem  4.1,  we  shall  postpone  the  proof  of  this 
theorem,  but  we  shall  now  show  its  relationship  to  reproducing  mdf's,  for 
which  the  following  lemma  will  be  needed. 

Lemma  5.3.  Let  functions  K  and  H  satisfy  the  same  hypotheses  as  in 
Lemma  4.2  except  that  (i)  K  vanishes  at  the  integers,  and  (ii)  6r+1H  is 
given  in  (o+p,  0-p)  by  a  spline  function  of  degree  2m  -  1  with  knots 
as  specified  in  condition  (iii)  of  Theorem  5.2.  Then  (4.2)  holds. 

Proof .  The  proof  is  the  same  as  that  of  Lemma  4.2  down  to  equation  (4.3) 
except  that  q  is  a  spline  of  degree  m  -  1  with  knots  as  specified  in 
condition  (ii),  and  it  follows  from  the  expression  (3.5)  for  G  that  it 
vanishes  at  the  knots  of  q.  Thus  m-fold  integration  of  (4.3)  gives 
O  -  (-1)m+r+1E  G(x)[q(m‘1)(x  +  0)  -  q(m_1)(x  -  0)1  , 


where  the  summation  is  over  the  knots  of  g  in  |a  +  p,  S  •  p).  Since  G 
vanishes  at  these  knots,  c  »  0,  as  required.  □ 

Theorem  5.4.  The  unique  interpolation  formula  determined  by  Theorem  5.2 
is  the  unique  mdf  of  the  class  ^xrm* 

Proof .  The  proof  is  identical  to  that  of  Theorem  4.3,  except  that  the 
role  of  Lemma  4.2  there  is  now  assumed  by  Lemma  5.3.  □ 


6.  COMPACT  EXPRESSION  FOR  mdf  BASIC  FUNCTIONS .  If  there  is  a  formula 
of  Fj  with  a  basic  function  that  satisfies  conditions  (i)  and  (ii)  of 
Theorem  4.1,  this  basic  function  can  be  regarded  as  a  spline  function  of 
degree  d  with  multiple  knots  of  multiplicity  d  -  m  +  1 .  In  general , 
therefore,  it  has  a  unique  expression  (see  [2])  of  the  form 

d  n  .  . 

L(x)  -  I  £  [c  (x  “  °  +  -  0  +  j)+]  ,  (6.1) 

i=m  j»0  3  3 

where  n  denotes  the  largest  integer  contained  in  V  and  y^  -  max(y^,0). 
The  coefficients  c^  and  g^  are  subject  to  the  constraints  arising  from 
Theorem  3.1,  which  can  be  wr itten  as 


(a)  I  j^,  =  0 
j-0  J 


at 

(b)  I  jkgA1  -  0 
j-0 


(k  -  0,1,  ...,r;  i  «  m,m  +  1,...,d) 


(6.2) 


If  W  is  an  integer,  the  second  term  of  the  summand  in  (6.1)  is  absent  (as 
are,  of  course,  the  constraints  (6.2)(b))« 

For  a  formula  of  t5iat  satisfies  conditions  (i)  and  (ii)  of  Theorem 

5.2,  in  general  there  must  be  added  to  the  right  member  of  (6.1)  the 


expression 


I  V*  -  j)?"1  ,  (6.3) 

jez 


where  E  denotes  the  set  of  integers  contained  in  I,  and  the  constraints 

I  jkh..  -  0  (k  -  0,1,  ...,r)  (6.4) 

j8E 


must  be  satisfied.  However,  when  a  or  B  is  an  integer,  the  addition  of 
(6.3)  is  not  required. 

We  shall  now  show  that  by  taking  into  account  condition  (iii)  of  Theorems 
4.1  and  5.2  and  by  introducing  certain  special  spline  functions,  we  can 
rewrite  (6.1)  in  a  form  involving  a  much  smaller  number  of  undetermined 
coefficients  and  can  also  avoid  the  necessity  of  considering  separately  the 
cases  in  which  P,  a,  or  8  is  an  integer.  For  this  purpose  we  shall  need 
the  following  lemma. 

Lemma  6.1 .  For  every  nonnegative  integer  r  and  every  positive  integer 
n  ^  r  +  1 ,  there  is  a  unique  polynomial  prn  &  such  that 


j-1 


JkPr„')> 


Ok 


(k  -  0,1, 


,r) 


(6.5) 


Proof .  This  follows  easily  from  well  known  properties  of  orthogonal 
polynomials,  but  it  is  also  readily  seen  as  follows.  Equations  (6.5)  may  be 
regarded  as  a  system  of  r  +  1  linear  equations  in  the  r  +  1  coefficients 
of  Prn«  The  latter  system  has  a  unique  solution  if  and  only  if  the  corre¬ 
sponding  homogeneous  system  has  only  the  trivial  solution.  But  any  solution 
of  the  homogeneous  system  gives  rise  to  a  polynomial  P  e  if  such  that 

n 

l  P(j)q(j)  -  0  (6.6) 

3-1 
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for  all  q  6  In  particular,  one  may  take  q  ■  P,  so  that  (6.6)  beoomes  a 

sum  of  squares,  and  therefore  P  vanishes  for  j  »  1,2,...,n.  Since  n  >  r, 

P  is  identically  zero .  □ 

We  now  define  the  splines  of  degree  r, 

n 

srn(x)  =  x+  -  I  Prn(3)(x  “  * 

j-1 

n 

S?  (X)  =  x£  -  l  Prn ( j ) ( x  ♦  j)  +  . 

j-1 

Parenthetically,  we  remark  that  by  means  of  (6.5)  and  the  identity 

y+  -  yr  -  (-i)r(-y)+ 

it  is  easily  shown  that  SJn(x)  -  (-1 )r+1Srn(-x)  .  We  observe  also  that  Srn 
has  its  support  in  (0,n)  and  Sjn  in  (-n,0)  . 

Note  that  condition  (iii)  of  Theorem  4.1  or  5.2  implies,  in  general,  that 
certain  knots  that  the  function  or+1L  would  otherwise  be  expected  to  have 
are  absent  (or  reduced  in  multiplicity  in  special  cases  of  a  reproducing 
formula).  Using  the  notation  of  (6.1),  this  means  that 

Ar+1Cij  *  Ar+1g^  =  0  <j  =  1  #2, . .  .,N  -  r  -  1>  i  *  m,m  +  1,...,d)  , 
where  the  finite  differences  are  taken  with  respect  to  j,  and  N  is  the 
largest  integer  less  than  V.  (Note  that  N  differs  from  n  when  W  is  an 
integer.)  This  implies  the  existence,  for  i  *  m,m  +  1,...,d,  of  a 
polynomial  such  that  c^  =  qi < j )  for  j  -  1,2,...,N  but  not,  in 

general,  for  j  =  0.  Similar  remarks  apply  to  9^ j •  We  conclude  from  these 
facts  and  the  constraints  (6.2)  that 
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(6.7) 


l  k.fx  -  «  -  j>+  <3tJ*  -  8  ♦  j>*] 

j-0  J  J 

*  cisr5”i,(x  “  0)  +  9i8rNr”i)(x  " 

(1  *  m,m  +  1 i • * d)  / 

where  c^  and  are  obtained  by  multiplying  c^g  and  g^Q  by  appropriate 

constants  depending  only  on  i  and  r ,  and,  in  case  i  >  r,  a  derivative  of 
negative  order  denotes  that  particular  integral  of  corresponding  order  that 
varnishes  identically  for  x  <  a.  A  little  reflection  will  convince  the 
reader  that  the  right  member  of  (6.7)  is  a  valid  substitution  even  when  U  is 
an  integer  (and  N  *  U  -  1).  However,  in  the  case  of  a  reproducing  formula, 
the  expression  (6.3)  now  must  always  be  added,  even  though  a  or  6  (or 
both)  is  an  integer . 

Accordingly,  (6.1)  can  be  expressed  uniquely  in  the  form 

L(x)  -  f  [c^J^U  -  a)  +  g18^r"i>(x  -  B)]  .  (6.8) 

i-m 

A  function  L  satisfying  conditions  (i)-(iii)  of  Theorem  4.1  that  is  the 
basic  function  of  a  formula  of  FIrm  has,  therefore,  a  unique  expression  of 
the  form  (6.8).  Similarly,  a  function  L  satisfying  oonditions  (i)-(iii)  of 
Theorem  5.2  that  is  the  basic  function  of  a  formula  of  has  a  unique 

expression  of  the  form  (6.8)  with  (6.3)  added. 

The  coefficients  c^  and  g^  (and  hj  in  the  reproducing  case)  must 
satisfy  certain  constraints.  These  will  now  be  described. 

(a)  In  the  reproducing  case  the  r  +  1  conditions  (6.4)  must  be 
satisfied.  When  this  is  the  case,  the  expression  (6.3)  vanishes  for  x  £  0. 

(b)  The  function  given  by  (6.8)  has  its  support  in  I  when  d  -  r. 

This  is  also  true  in  the  reproducing  case  if  the  preceding  condition  (a)  is 
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fulfilled.  However,  when  d  -  2m  -  1  >  r,  (6.8)  gives,  for  x  >  6,  e 
polynomial  of  degree  d  -  r  -  1 .  The  vanishing  of  the  d  -  r  coefficients  of 
this  polynomial  involves  d  -  r  constraints. 

(c)  in  the  general  case,  exactness  for  the  degree  r  requires  that 
the  r  +  1  conditions  (5.1)  be  satisfied. 

(d)  In  the  reproducing  case,  (c)  is  replaced  by  conditions  (1.2).  This 
involves  t  effective  constraints. 

(e)  The  form  of  (6.8)  ensures  that  6r+1L  shall  be  given  in 

(a  -f  p,  6  -  p)  by  a  polynomial  in  the  general  case,  and  by  a  polynomial 
except  for  discontinuities  in  the  (2m  -  1)th  derivative  in  the  reproducing 
case.  However,  in  the  case  when  d  -  r  >  2m  -  1 ,  condition  (iii)  of  Theorem 
4.1  or  5.2  involves  also  a  reduction  in  degree  from  d  to  2m  -  1 .  This 
constitutes  d  -  2m  +  1  constraints. 

7.  PROOFS  OF  THE  CHARACTERIZATION  THEOREMS.  We  now  have  the  machinery 
needed  to  prove  Theorems  4.1  and  5.2. 

Proof  of  Theorem  4.1.  If  there  is  a  formula  of  FjrlB  whose  basic 
function  satisfies  conditions  (i)-(iii)  of  Theorem  4.1,  that  function  has  a 
unique  expression  of  the  form  (6.8)  with  parameters  c^  and  g^  satisfying 
conditions  (b),  (c),  and  (e)  of  the  preceding  section,  to  the  extent  these 
conditions  are  applicable.  On  the  other  hand,  if  there  is  an  expression  of 
this  form  with  parameters  satisfying  these  conditions,  then  it  is,  in  fact, 
the  basic  function  of  such  a  formula . 

Now,  (6.8)  contains  2(d  -  m  +  1)  undetermined  parameters.  The  numbers 
of  constraints  involved  in  conditions  (b),  (c)  and  (e)  are,  respectively, 
d  -  r,  r  +  1  and  d  -  2m  +  1 .  (Note  that  the  integers  d  -  r  and 
d  -  2m  +  1  are  not  both  different  from  zero.)  The  total  number  of 


constraints  is  2(d  -  m  +  1),  ths  sans  as  ths  number  of  parameters  •  Without 
spelling  out  the  constraints  in  detail ,  it  is  easily  verified  that  they  are 
linear  equations  in  the  parameters •  The  parameters  must  therefore  satiefy 
2(d  -  m  +  1)  equations  in  as  many  unknowns*  lb  prove  Theorem  4*1  it  is 
sufficient  to  show  that  this  linear  system  is  nonsingular . 

This  is  the  case  if  the  corresponding  homogeneous  system  has  only  the 
trivial  solution.  In  fact,  the  only  one  of  the  constraint  equations  that  has 
a  number  other  than  0  on  its  right-hand  side  is  the  one  obtained  by  taking 
i  *  0  in  (5.1).  Thus,  a  function  K  of  the  form  (6.8)  whose  parameters 
satisfy  the  homogeneous  system  has  the  property  that  TR  annihilates  *  . 

This  function  K  therefore  fulfills  the  requirements  for  both  K  and  H  in 
Lemma  4.2.  Consequently,  by  Lemma  4.2, 

m 

!  lK(B,(x)J2dx  -  0  . 

mm 

By  the  same  reasoning  used  in  the  proof  of  Theorem  4.3,  it  follows  that  K  Is 
identically  zero.  Thus,  the  homogeneous  system  has  only  the  trivial  solution. 

Proof  of  Theorem  5.2.  The  perceptive  reader  may  have  noticed  that  the 
possibility  of  m  ■  0,  though  allowed  in  Theorem  4.1,  is  excluded  in  Theorem 
5.2.  In  fact,  a  reproducing  mdf  with  m  -  0  is  somewhat  meaningless,  for  the 
following  reason.  Application  to  this  case  of  the  criteria  that  we  have 
developed  would  lead  to  a  solution  in  which  the  basic  function  of  the 
corresponding  mdf  without  the  reproducing  requirement  is  modified  by 
arbitrarily  assigning  at  the  integers  the  values  given  by  (1.2),  even  though 
these  are  inconsistent  with  the  values  at  neighboring  arguments.  Thus,  the 
resulting  basic  function  would  have  removable  discontinuities  at  the  Integers 


in  I.  Strictly  speaking,  such  a  function  is  not  piecewise  continuous,  and 
therefore  is  not  the  basic  function  of  a  formula  of  the  class  FIr0. 

If  there  is  a  formula  of  Fj*£  whose  basic  function  satisfies  aonditions 
(i>— (iiij  of  Theorem  5.2,  that  function  has  a  unique  expression  of  the  form 
(6.8)  with  (6.3)  added,  and  the  parameters  c^,  ,  and  hj  satisfy 

conditions  (a),  (b),  (d),  and  (e)  of  the  preceding  section,  to  the  extent 
these  conditions  are  applicable.  On  the  other  hand,  if  there  is  an  expression 
of  this  form,  with  parameters  satisfying  these  conditions,  then  it  is,  in 
fact,  the  basic  function  of  such  a  formula. 

Now  (6.8)  and  (6.3)  together  contain  2(d  -  a  +  1)  ♦  t  undetermined 
parameters  •  The  number  of  constraints  involved  in  aonditions  (a),  (b) ,  (d), 
and  (e)  are  respectively,  r  ♦  1,  d-r,  t,  and  d  -  2m  +  1.  The  total 
number  of  constraints  is  2(d  -  m  +  t)  +  t,  the  same  as  the  number  of 
parameters.  As  in  the  general  case,  all  the  constraints  are  linear  equations 
in  the  parameters,  and  they  constitute  a  linear  system  having  a  square 
coefficient  matrix. 

The  remainder  of  the  proof  is  the  same  as  for  Theorem  4.1,  except  that 
Lemma  5.3  now  assumes  the  role  played  by  Lemma  4.2  in  the  earlier  proof. 

8.  SOME  mdf's  ARE  PREVIOUSLY  PUBLISHED  FORMULAS.  In  some  instances  the 
min imized-der ivatlve  formula  of  a  class  turns  out  to  be  a  previously  published 
formula.  Table  I  lists,  for  the  cases  known  to  us,  the  class  Flrjn  or 
Fl^m  involved,  the  name  of  the  originator,  the  publication  citation,  and  the 
year  of  publication .  Two  of  the  papers  cited  contain  a  large  number  of 
formulas,  and  in  these  cases  the  particular  formula  is  identified.  In  two 
instances  in  which  the  published  formula  contains  an  unspecified  parameter, 
the  numerical  value  of  the  parameter  that  yields  the  mdf  is  given  in  a 
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footnote .  The  entry  "Both"  in  the  fourth  column  means  that  FIrm  and 
are  identical  for  the  case  involved. 


TABLE  I.  Previously  Published  Formulas  that  are  mdf's 


I 

r 

m 

Rep  or 
Nonrep 

Originator 
and  Citation 

Publicatior 

Year 

(-2,  2) 

1 

2 

Nonrep 

Jenkins  [5] 

1927 

(-2  ,  2) 

2 

2 

Both 

Xarup  [6] 

1898 

(-2,  2) 

2 

3 

Both 

Greville  [1]  (105) 

1944 

(-5/2 ,  5/2) 

2 

2 

Nonrep 

Greville  [1]  (67 )2 

1944 

(-5/2,  5/2) 

3 

2 

Nonrep 

Greville  (1)  (69) 

1944 

(-3,  3) 

3 

2 

Nonrep 

Greville  (1]  (73)3 

1944 

(-3,  3) 

3 

3 

Nonrep 

Vaughan  [10]  "C" 

1946 

(-3,  3) 

3 

2 

Rep 

Henderson  [4] 

1906 

(-3,  3) 

4 

2 

Both 

Shove lton  [8] 

1913 

(-3,  3) 

4 

3 

Both 

Sprague  [  9] 

1880 

2 With  a13  =  13/80. 

3With  aQ4  -  -7/108. 
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